Josephson effect through an anisotropic magnetic molecule 
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We study the Josephson effect through a magnetic molecule with anisotropic properties. Perform- 
ing calculations in the tunneling regime, we show that the exchange coupling between the electron 
spin on the molecule and the molecular spin can trigger a transition from the tt state to the state, 
and we study how the spin anisotropy affects this transition. We show that the behavior of the 
critical current as a function of an external magnetic field can give access to valuable information 
about the spin anisotropy of the molecule. 

PACS numbers: 74.78.Na, 74.45. +c, 73.63. Kv, 75.20.Hr 



I. INTRODUCTION 

The Josephson effect 1 -^ is a striking manifestation of 
many body physics and macroscopic quantum coherence 
in condensed matter systems. While early investigations 
concerned mainly bulk superconducting junctions sepa- 
rated by an insulating barrier, in the last decades it has 
become a very active field of study in the context of meso- 
scopic physics. Indeed the insulating barrier can be re- 
placed by a conductor or a nano-device that can be as 
small as a quantum dot or a single molecule. In this 
sense the study of the Josephson current can provide a 
novel way to investigate the electronic properties of the 
nano-object, which is sandwiched between the supercon- 
ducting electrodes. More than a decade ago, it was pre- 
dicted using the Krein theorem 3 -"— that when a singly 
occupied quantum dot in the Coulomb blockade regime 
is inserted between the superconductors, the Josephson 
current phase relation acquires a 7r shift, i.e., the critical 
current has the opposite sign from that of a tunnel junc- 
tion. A phase diagram of the tt-0 phase transition was de- 
rived later on for contacts with arbitrary transparency us- 
ing a combination of Hubbard-Stratonovich and saddle- 
point approximation. 6 Experimentally (for nanoscale de- 
vices) it was measured in superconductor-nanotube- 
superconductor systems.— This picture gets more com- 
plicated when the Kondo temperature is lower than the 
superconducting gap: a 0-junction state is restored^ 
albeit with a different current phase relationship. 

In recent years theoretical and experimental studies 
have addressed transport geometries where a molecule — 
artificial or otherwise - is inserted between two 
electrodes.— ~— This goes one step beyond the study of 
transport through quantum dots because the molecule 
has internal degrees of freedom (such as vibrations and 
possibly spin). On the one hand, such degrees of free- 
dom have an effect on the electronic current, on the other 
hand, the current itself can be considered as a probe of 



the inherent mechanisms of the molecule. 

A subfield of molecular electronics is called molec- 
ular spintronics: it focuses on molecules which have 
an intrinsic spin ) 16 ' 17 and it is expected that electron 
transfer through the molecule can trigger changes in the 
molecule spin because of the existence of an exchange 
coupling with the electron spin. Such molecules (such 
as a buckminsterfullerene doped with a magnetic atom) 
may have an isotropic spin, or otherwise the spin may 
have a preferred direction due to the crystalline struc- 
ture of the molecule (this is the case of Mni2 acetate) . 
Recently, there have been some efforts to describe and 
measure transport through molecular spintronics devices 
with normal metal or ferromagnetic leads,— ~— with an 
emphasis on master equations approach on the theoret- 
ical side. Nevertheless, efforts in the field of molecular 
spintronics with superconducting electrodes are still at 
their beginning stage. 

A recent theoretical work focused on the Josephson 
current through an isotropic magnetic molecule, via per- 
turbative calculations in the tunneling Hamiltonian as 
well as numerical renormalization group calculations. 22 
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FIG. 1: (Color online) A magnetic molecule (e.g., M@C80, 
Mni2, . . .) connecting two superconductors via tunnel barriers 
£ L and i R . The exchange coupling between the molecular spin 
and the electronic spin can strongly modify the Josephson 
current. 
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It allowed us to draw a complete phase diagram of the 
7T-0 phase transition. An equivalent study of supercur- 
rcnt through molecules which have an anisotropic spin, 
which magnetization can tunnel, and which are subject 
to a (weak) external magnetic field is still lacking. This 
is precisely the focus of the present work. One of the 
challenges of this work is that we have to deal with a 
large number of parameters: the exchange coupling J 
between the dot electron spin and the molecule spin, the 
anisotropy constant D, and the coefficient B 2 for quan- 
tum tunneling of magnetization of the molecule, the de- 
pendence on external magnetic field B as well as the dot 
level £d, which can be adjusted by a gate voltage. Note 
that it is now experimentally possible to manipulate the 
anisotropy parameters of magnetic molecules . 23 ' 24 One 
of our goals is to determine to what extent the measure- 
ment of the critical current can provide information of the 
sign or magnitude of such parameters. For simplicity, we 
focus on the regime where the superconducting gap is 
much larger than the Kondo temperature, which allows 
us to focus on weak coupling (small tunneling Hamilto- 
nian) calculations. Also, we restrict the analysis on the 
simplest case of a molecule spin S = 1 to demonstrate 
the effect where the two main contributions due to spin 
anisotropy (easy axis anisotropy and quantum tunneling 
of magnetization) are present. 

The outline of the paper is as follows. In Sec. [Ill we in- 
troduce the model for the magnetic molecule connected 
to two superconducting leads and we compute the ex- 
pression of the Josephson current through this molecule. 
In Sec. IIII1 we study the effect of the anisotropic param- 
eters and of the adjustable experimental parameters on 
the sign of the critical current. Finally, we conclude in 
Sec. HV] 



II. MODEL 

A. Hamiltonian 

The total Hamiltonian of the system (see Fig. [1} con- 
sists of the three terms H = + T-L s +1-Lf The first one 
is the Hamiltonian of the molecule, 
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Hd = H m + e d dU a + Ud\d^d\d v 



(1) 



where ed is the electronic level of the molecule implied in 
the transport, and U is the Coulomb interaction strength. 
The 4 and d a are electronic creation and annihilation 
operators on the electronic level in the molecule. Since U 
is typically much larger than the other energies in the sys- 
tem, we consider the limit of infinite Coulomb interaction 
U — > oo, thus only one electron is allowed to occupy the 
dot. With this assumption, the Hamiltonian % m , which 
characterizes the magnetic properties of the dot, reads 
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FIG. 2: (Color online) Spin states for S — 1 spin, and cou- 
pling between these states due to the different terms of the 
Hamiltonian. B2 induces tunneling between |l) m and j — l) m 
states; J and B induce |l) m o |0) m and |0) m o | — l) m 
tunnelings. 



where S z is the molecular spin and s z the spin of the 
electron on the molecule (if present). J is the exchange 
coupling between molecular and electronic spin, D > 
is the easy axis anisotropy constant, B2 is the coeffi- 
cient of quantum tunneling of magnetization (QTM), 
and B is the external magnetic field. Figure [2] shows 
how these terms couple the states of the molecule in 
the case of a spin 5=1. In order to avoid a too 
large number of parameters, we have made some sim- 
plifying assumptions when writing this Hamiltonian: the 
anisotropy terms are not affected by the charge of the 
dot level (this should be the case for systems like M@Cgo, 
but not for molecules like Mni2)r^ the magnetic held is 
taken parallel to the spin anisotropy,^, and higher or- 
der terms (— _B 2 n/2)(S'+™ + S 2 ? 1 ) are neglected (they are 
usually small). 

The second term corresponds to the superconducting 
parts, described by the BCS Hamiltonian 



l,k,a 



t.k 



(3) 

where = h 2 k 2 /2m — E F is the dispersion relation 
for free electrons, c] fe and c £ k are electronic cre- 
ation and annihilation operators in the superconductors, 
I enumerates left (£ = L) and right (£ = R) leads, 
A L( R),fe = Ae ±lv ^ 2 , with A the superconducting gap 
and if the superconducting phase difference along the 
junction. 

The last term is the tunnel Hamiltonian between the 
leads and the molecule 



H t = \^^d\c l k a + H.c. 



(4) 
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where t^k ar e the tunneling amplitudes. By perform- 
ing a gauge transformation for t h> k = t^^e^^ 4 , t R ^ — 
t^,ke-~ lip ^ and simultaneously for c^^.a = c L ,fc,<xe~ Jip / 4 , 
CR,fc,(x = c Ri fe i<T e Jip / 4 , one can "move" the dependence on 
<p from Ai.k to t^k and C£.k,a in Eqs. ([3]) and We also 
perform a Bogoliubov transformation 2 to diagonalize the 
BCS Hamiltonian, which takes the following form: 



(S 2 + + S 2 _) + JSs, (2) 



n s = 
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E kl{k,all,k, 



(5) 



3 



and the tunneling Hamiltonian reads 



H t = 



[kkdl (u k \ k r7 + cr «fe7i,fe,-<r) + H - C - > ( 6 ) 



,k,a 



where 7I k and ^ye.k.a are the quasiparticle creation and 
annihilation operators, = + £k/Ek)/2 and = 
— ek/Ek)/2 are the electron and hole coefficients, 
and Ek — \J e 2 . + A 2 is the energy dispersion. In the 
following calculations we will consider for simplicity the 
case of symmetric contacts, thus t^^ — t R _k = tk- 



B. Specific Hamiltonian for the 5* = 1 case 

In the following, we will for simplicity restrict our cal- 
culations to the case of a molecular spin with 5 = 1, 
which is the smallest value where easy axis anisotropy 
(parameter D) and QTM (parameter B 2 ) are nontrivial. 
As the electron occupation of the level is restricted to 
or 1, let us write explicitly the molecule Hamiltonian in 
each case [see Eqs. (HJ and ©]. 

For the empty electronic level, we have Hd = H-m,o, 
and we use the basis {|0) o |l) m , |0) c |0) m , |0) c | - l) m }, 
where |0} e represents the empty electronic state and 
\S z ) m the states of the molecule with spin projections 
S z — 1. 0, —1. The matrix elements of H m are 



H n 



B-D -B 2 

-B 2 -B-D 



(7) 



The eigenvalues are noted Eq^ [i = 1,2,3), and the cor- 
responding eigenvectors are b;. Below we will use the 
matrix b^ = [b 1; b 2 , b 3 ], which consists of columns of 
eigenvectors (first index enumerates columns, the second 
enumerates rows), and the inverse matrix bij — (bij) . 

When the electronic level is occupied by one electron, 
we have Hd.i = H m ,i + e d, and we use the uncoupled 
spin basis |s) e |S , ;; ) m (with s =t,i and S z = +1,0,-1). 
The matrix representation of 7^1,1,1 can then be decom- 
posed as two independent 3x3 submatricess^i H m .i = 
diagfH+ ^"H" -J, with 



3B/2 + J/2- D -B 2 
-B/2 J/y/2 

-B 2 J/V2 -B/2- J/2 



D 



in the basis {|t)o|l)m, |l)c|0) m , |1%| - l)m} and 



(8) 



B/2- J/2- D J/y/2 -B 2 
'HZ i J/V2 B/2 

-B 2 -35/2 + J/2- D 

(9) 

in the basis {| |) c |l) m , | t)o|0) m , | |) c | - l) m }. These 
matrices have eigenvalues E^ { , and corresponding 

As previously, we define matrices 



eigenvectors a- , a, 

'af, a^, a^], and inverse matrices af- — (a*) 



C. Josephson current 

The Josephson current through the molecule can be 
calculated using perturbation theory in the tunneling 
Hamiltonian 'Ht'r the first nonvanishing term is given 
by 



I = ^--^(gs\H t (E gs - Uo)- x Me & 



x n t (E g 



Ho)- 1 
Hor'Htlgs), (10) 



where Hq = Hd + H s . The ground state |gs) is the 
occupied state with lowest energy, thus it has energy 
E gs = minji?^}, and |gs) = |a?), where i = 1,2,3 spec- 
ifies the state number and £ = ± is the block index. 
Note that the dot-lead coupling induces energy shifts for 
the occupied states of the dot, starting at order 2 in Ht- 
However, we do not need to compute these shifts, as they 
will be identical for the two single occupied states, and 
they can be included in the value of ed (see Ref. |H for a 
multilevel case where these shifts have to be computed). 

As was shown in Ref. y, in the absence of coupling 
to a molecular spin, the perturbative approach allows us 
to understand the ir state due to large Coulomb interac- 
tion on the dot: the order of the electrons of a Cooper 
pair is necessarily reversed during tunneling through the 
dot, which gives opposite sign for the current due to the 
singlet nature of the Cooper pair. Here, the exchange 
coupling between the electron spin and a molecular spin 
means that the occupied state of the dot is a linear com- 
bination of states involving in general both | f) and 1 1) 
states of the electron spin. This creates the possibility 
of spin-flip processes: a spin-up electron tunneling in the 
dot can tunnel out as a spin-down electron for example. 
With such a spin flip, it is now possible for a Cooper 
pair to tunnel through the dot without reversing the or- 
der of electrons, thus contributing to positive current. 
In the presence of exchange coupling with a molecular 
spin, one can thus expect that, among all the lowest- 
order processes contributing to the Josephson current, 
some of them will contribute to negative current, and 
some others to positive current. The global sign of the 
current will thus depend on the relative weight of the dif- 
ferent processes, which are a function of the parameters 
of the molecule Hamiltonian. 

Expressing in Eq. (fTO)) the action of the tunneling 
Hamiltonian on the eigenstates introduced in the pre- 
vious section, a lengthy but straightforward calculation 
gives eventually 
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FIG. 3: The dependence of the normalized critical current i c 
as a function of the exchange coupling J for an isotropic mag- 
netic molecule (D = B2 = 0) and dot level ed/A = —5 in the 
absence of magnetic field (B = 0). 



III. RESULTS AND DISCUSSION 

For reference, we start by analyzing Eq. (fTTj) as a 
function of exchange coupling J, when no anisotropy 
is present (D = B 2 = 0) and without magnetic field 
(B = 0). As shown in Fig.[3]Ja), the current is suppressed 
both by negative and positive J. For negative J (ferro- 
magnetic coupling) the system always remains in the n 
state (i c < 0). For positive J (antiferromagnetic cou- 
pling) a 7T-0 transition occurs for J/A ~ 10 (the precise 
value is slowly varying with e d ). This behavior can be 
understood by looking at the formula for the current^ 



Here 



(a g s,i & n + ags^siX&iiafj + ^13^-) 



Ek + Eqa — Egs — e d 



± (a gs ,ibi3 + a gs ,3&33)(b3ia^ + &33a^) 



^fe + ^0,3 _ E gs — e d 



Sfe + E 02 — Er, s — £d ' 



(13) 



where 



, and C = -C- Equations (jTTj) - (fT3"j) 

represent the main results of this paper. Because we have 
performed a lowest-order tunneling calculation, we get a 
simple I — I c shop dependence of the current. However, 
the study of value of the critical current I c (in addition 
to its sign) will give us precious information on the sys- 
tem. At zero temperature the sums over k and k! should 
be taken over the energy region ek,£k' > 0. Both sum- 
mations over k can be replaced by the integration over 
energy e: J2k I d e P( € )> where p(e) is a density of 
states. 

Our formulas of course contain the known result for the 
case where there is no molecular spin and no magnetic 
field (B — B2 = J = D = 0)>22 we obtain a 7r-junction 
with negative critical current, 



7(0) 
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h f£ (E k + E k '){E k - e d )(E k , - e d ) 



4e r L r R A 2 

h 4tt 2 J E 1 E 2 {Ei+E 2 ){E 1 -e A ){E 2 -e A )' 


(14) 



where we assume constant density of states p(e) 



Po 



2m/ nh , tunneling rates T 



L(R) 



KPotl iRV and E 1{2) = 



A 2 



-~1(2)- 



In the next section we analyze the dependence of the 
dimensionless critical current i c = I c /\Ic°^\ on the mag- 
netic molecule parameters J, D, B 2 , dot energy e d , and 
external magnetic field B. Positive i c > corresponds to 
the 0-junction phase, negative i c < corresponds to the 
7r-j unction phase. 
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3£fc£fc' [3J/2 + Ek + Ek> Ek + Ek 



(15) 



where £k = Ek + J — e d - The first term depicts the 
transfer of a Cooper pair involving a change of the total 
coupled spin (electronic and molecule) during the inter- 
mediate state [e.g., see Fig. EJa)], while the second term 
corresponds to a Cooper pair without change of total spin 
during the intermediate state [e.g., see Fig. 0Jb)]. For 
large positive J, the first term becomes smaller than the 
second one, and the sign of the current changes, which 
explains the tt-0 transition. 

Note that there is no change of ground state associated 
with this transition occurring for large positive J, hence 
the critical current shows a smooth change from negative 
to positive value, passing continuously through arbitrary 
small values. This is to be contrasted with 0-tt transition, 
which is due to the crossing of energy levels leading to 
a change of ground state ^ 30 ' 31 where an abrupt change 
of the critical current can be observed (see, e.g., Figs. 4 
and 6 in Ref.[H). 

We will now consider the effect of the anisotropy (D 
and B 2 ) and of the magnetic field B on the critical cur- 
rent, especially near the 7T-0 transition. We assume that 
the superconducting gap is independent of the magnetic 
field. Figure [3](a) shows the effect of D and B 2 on the 
transition; the surface shows the values of the parameter 
for which the current is zero. Above the surface the sys- 
tem is in the 7r-junction phase (i c < 0), while under the 
surface the system is in the zero phase (i c > 0). One can 
see that both D and B 2 move the 7T-0 transition to higher 
values of J. This is confirmed by Figs. [5jb) andOJc), 
which correspond to cuts of the three-dimensional (3D) 
plot for fixed values of B 2 and D, respectively. On these 
panels, the different curves correspond to different val- 
ues of the magnetic field B: we see that increasing the 
magnetic field tends to push the system toward the 0- 
j unction phase (note that the results are insensitive to 
the sign of B). On the 3D plot Fig. [5ja) the effect of the 
magnetic field B is thus to shift the zero current surface 
as shown with magenta arrows, and also to somewhat 
smear the sharp behavior in B 2 as shown with blue ar- 
rows. 
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FIG. 4: (Color online) Illustration of two typical tunnel- 
ing processes leading to the transfer of a Cooper pair. The 
presence of strong Coulomb interaction prohibits the double 
occupation on the dot and the electrons are transferred one 
by one. Because of the exchange coupling J between the elec- 
tron spin and the molecular spin, the state of the occupied 
dot (black circle on the figure) is characterized by the total 
spin, S — 1/2 and S + 1/2. These two levels are separated 
by an energy 3J/2. The process where the intermediate state 
of the occupied dot (a) is different from the initial one and 
the process where the intermediate state is the same as the 
initial one (b) contribute with different signs to the Josephson 
current. The competition between these two processes leads 
to the existence of the 7T-0 transition; see Eq. (|15[) . 



Up to now, we have studied the phase diagram of the 
system as a function of the exchange coupling J and of 
the anisotropy parameters D and i?2- However, for a 
given molecule, these parameters have usually a fixed 
value. We will now study the behavior of the critical 
current when the experimentally adjustable quantities, 
the external magnetic field B and the dot level ed, are 
varied. The goal is to understand how the values of the 
exchange coupling and of the anisotropy parameters will 
modify the behavior of the current as a function of B 
and £d • This could be an original way to obtain informa- 
tion on the exchange coupling and on the spin anisotropy 
in the molecule, by measuring the critical current of the 
tunnel junction and varying B and ea- 

The different panels of Fig. [5] show the behavior of the 
critical current as a function of B and for various values 
of the dot level e<j, the exchange coupling J and of the 
anisotropy parameter D (for simplicity, we have taken 
B2 =0). Each column is for a given value of J: deep in 
the 7r-junction regime J / A = —20 (left), in the intermedi- 
ate regime J/ A = 11 (middle), and deep in the 0-junction 
regime J/A = 20 (right). The top panel of each column 
is for D = 0, while the two bottom panels of each col- 
umn are for nonzero values of D as indicated. The richer 
behavior is obtained when the exchange coupling has a 
value that allows us to observe the 7T-0 transition, here 
in the second column for J/A = 11. Without anisotropy 
[Fig.[BJd)], we see that by sweeping the magnetic field we 



71-junction 




FIG. 5: (Color online) (a) n- and 0-junction regions as a 
function of the J, D, and -82- The surface divides 3D space 
(J,D,B2) to the top region in a 7r-junction regime and to 
the bottom one with 0-junction regime; at the surface current 
is zero. Magnetic field is zero, B — 0, its increasing leads 
to "shift" and "smear" of the surface as shown by magenta 
and blue arrows, (b) 7r-0 transition diagram in (J,D) space 
at B2 — 0. Different curves correspond to the different B's: 
B/A = 0.0 (solid), 0.5 (dashed), and 1.0 (dotted), (c) tt-0 
transition diagram in ( J, B2) space at different B and D = 0. 
The nonzero D and B2 increases the critical J c (see Fig. |3J) 
and magnetic field mainly decreases J c . 



can observe the 7T-0 transition. In the presence of small 
anisotropy [Fig.|5Je)], we observe a nonmonotonic behav- 
ior as a function of B, with the modulus of the critical 
current \i c \ decreasing as a function of B for small B, 
but increasing for large B. Finally, for larger anisotropy 
[Fig-EJf)] j I j c | i s everywhere increasing as a function of B. 
Note that, between panel Fig. M,d) (D — 0) and panel 
Fig-EJf) (D/A = 4), the order of the curves as a function 
of £d has been reversed. When J is much larger than the 
superconducting gap (right column, with J/A = 20), the 
system is deep in the 0-junction phase, but the anisotropy 
has a visible impact on the curves: comparing Fig. [6jg) 
(for D/A = 0) with Figs. H^h) and Mi) (for D/A = 7 
and 14) , we see that when D is large enough, the slope of 
the critical current is the opposite of the one for small D. 
This is a consequence of the 7T-0 transition, which hap- 
pens for larger D. Finally, for negative J [Figs.|Ha)-[(){c) 
with J/A = —20], the anisotropy does not bring any 
qualitative change to the behavior of the current as a 
function of the magnetic field, and \i c \ always decrease 
with B. 

From the different curves shown in Fig. M we can de- 
duce that when J is positive (the antiferromagnetic cou- 
pling case), the anisotropy has a visible impact on the 
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FIG. 6: (Color online) Critical current as a function of the magnetic field B for different values of the exchange coupling J and 
the anisotropy parameter D (B2 = 0). The different curves in a single plot are for various dot levels ea: £d/A = —12.5 (solid 
cyan line), —10.0 (dashed green line), —7.5 (dotted red line), and —5 (dashed-dotted magenta line); all currents are normalized 
by the critical current obtained for A = —5 with B — J — D — B2 = [see Eq. (|14p ]. Left column: J/ A — —20; the system 
is deep in the 7r-junction regime, and the anisotropy parameter D does not change the curves qualitatively. Middle column: 
J/ A = 11; the system is near tt-0 transition, and the value of D has a great impact on the behavior of the curves: it can 
produce nonmonotonic behavior as a function of B [panel (e)], or reverse the slope of the curves compared to D = [panel 
(f)]. Right column: J/ A — 20; the system is in the 0-junction regime and D still has a visible impact, as it can change the 
slopes of the curves [panel (h)]. This is due to the presence of a (large) critical value D c above which the system is again in 
the 7r -junction phase (not shown). 



behavior of the critical current as a function of B, as 
it can produce a nonmonotonic behavior close to the tt— 
transition, and reverse the slope of \i c \ as a function 
of B when J is much larger than the critical value. On 
the other hand, for negative J (ferromagnetic coupling), 
the anisotropy does not have a qualitative effect on the 
critical current, and it merely reduces the value of |i c |. 



IV. CONCLUSION 

We have computed the Josephson current through a 
magnetic molecule in the tunneling regime, studying the 
effect of the exchange coupling with the molecular spin, 
and the spin anisotropy of the molecule. Performing 
a perturbative calculation starting from a Hamiltonian 
model, we have shown that an antiferromagnetic cou- 
pling between the electron spin and the molecular spin 
can induce a tt-0 transition. We have described how the 
spin anisotropy D and the quantum tunneling of magne- 
tization term Bi affect the transition. 

We have shown that by studying the behavior of the 



critical current as a function of the magnetic field and 
the level position (which are both experimentally tunable 
parameters), it is possible to get information on the value 
of the spin anisotropy D, even outside the range of the 
7T-0 transition. 

This work could be extended in several directions. The 
calculations could be performed for a larger molecular 
spin (albeit at the cost of heavier expressions) . One could 
also use anisotropy parameters that depend on the charge 
state of the molecule (and thus on the occupation of the 
dot in our model), which could describe more faithfully 
molecular magnets like Mni2*2£ One could also consider 
the case of an external magnetic field aligned along an 
arbitrary direction (and not along the anisotropy axis of 
the molecule), in order to describe experiments where 
it is not possible to control the anisotropy orientation. 
Such a magnetic field should have a strong impact on 
the current, as it will mix efficiently all the molecular 
states^ 

Finally, new possibilities could open up if one considers 
explicitly the Josephson current between type II super- 
conductors. In this case, it could be possible to con- 
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trol the value of the superconducting gap A with the 
applied magnetic field. Going to very small A would 
give large values of J/A, D/A, etc., and a very large 
parameter range of the system, including the 7r-0 transi- 
tion for J > 0, could be explored. In the same manner, 
it is possible to enhance the critical temperature T c and 
the second critical field H C 2 by decreasing the thickness 
of the superconductor^ This could allow us to use large 
values of the magnetic field. 
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